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Abstract: We present a general analysis of SU(3) breaking effects in the 
semi-inclusive tau hadronic width. The recent ALEPH measurements of the inclu- 
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sive Cabibbo-suppressed decay width of the r and several moments of its invariant 
mass distribution are used to determine the value of the strange quark mass. We 
obtain, in the MS scheme, m s (M^) = (119 ± 24) MeV to 0(al), which corresponds 
to m s (lGeV 2 ) = (164 ± 33) MeV, m s (4GeV 2 ) = (114 ± 23) MeV. 
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1. Introduction 



Among the free parameters of the Standard Model, the quark masses are the ones 
less precisely known. The lack of accurate measurements sensitive to quark mass 
effects and the theoretical uncertainties associated with the non-perturbative nature 
of QCD in the infrared region make quite difficult to perform reliable determinations 
of quark masses. In particular, the value of the strange quark mass has been a subject 
of great controversy in recent years. 

In the last version of the Review of Particle Physics [1], the running strange 
quark mass at 2 GeV in the MS scheme is quoted to be in between 60 MeV to 170 
MeV. This wide range reflects both the uncertainties in the hadronic input needed 
in strange quark mass determinations within the context of QCD Sum Rules and the 
spread in m s values obtained within lattice QCD calculations. 

The high precision data on tau decays [2] collected at LEP and CESR provide 
a very powerful tool to analyse strange quark mass effects in a cleaner environment. 
The QCD analysis of the inclusive tau decay width has already made possible [3] an 
accurate measurement of the strong coupling constant at the r mass scale, a s (M^), 
which complements and competes in accuracy with the high precision measurements 
of <Xg(M|) performed at LEP. More recently, detailed experimental studies of the 
Cabibbo-suppressed width of the r have started to become available [4, 5], which 
allows to initiate a systematic investigation of the corrections induced by the strange 
quark mass in the r decay width. First theoretical studies were presented in [6, 7]. 

What makes a m s determination from r data very interesting is that the hadronic 
input does not depend on any extra hypothesis; it is a purely experimental issue, 
which accuracy can be systematically improved. The major part of the uncertainty 
will eventually come from the theoretical side. However, owing to its inclusive char- 
acter, the total Cabibbo-suppressed tau decay width can be rigorously analyzed 
within QCD, using the Operator Product Expansion (OPE). Therefore, the theo- 
retical input is in principle under control and the associated uncertainties can be 
quantified. 

In the following we will compile and analyze in detail all what is presently known 
about quark mass corrections to the quark current correlation functions relevant in 
tau decay. In particular, we will investigate the size of these effects in the tau decay 
width and related observables, and the theoretical uncertainties of the corresponding 
predictions. 

Even with the relatively large uncertainties one gets from the present data, we 
will show that the strange quark mass determination from tau decays has already an 
accuracy good enough to reduce substantially the range quoted by the Particle Data 
Group. The foreseen r-charm and B-factories will further increase the precision of 
this measurement, allowing for much more detailed studies. Clearly, the tau decay 
data will provide in the future a precise determination of the strange quark mass 
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within QCD. 



2. Theoretical Framework 



The theoretical analysis of the inclusive hadronic tau decay width [8, 9, 10] involves 
the two-point correlation functions 



H% v (q) ee i J d*xe^ (0|T { V *(x) V^} |0> , 

0\T{A»(x)A^(0y}\0), 



i I d x e 



(2.1) 
(2.2) 



associated with the vector, Vfi{x) = q^qi, and axial-vector, A^(x) ee q^^qi, 
colour-singlet quark currents. The subscripts i,j denote the corresponding light 
quark flavours (up, down, and strange). These correlators have the Lorentz decom- 
positions 



LL ij,V/A 



(?) = (-<r q 2 + <t<f) nJj,v/A(q 2 ) + ng, VM (g 2 ) , (2.3) 



where the superscript in the transverse and longitudinal components denotes the 
corresponding angular momentum J = 1 (T) and J = (L) in the hadronic rest 
frame. 

The imaginary parts of the two-point functions H^ v , A (s) are proportional to 
the spectral functions for hadrons with the corresponding quantum numbers. The 
semi-hadronic decay rate of the r lepton, 



r r 



hadrons v T {l) 



e v e v T {i) 



(2.4) 



can be expressed as an integral of the spectral functions ImIl T (s) and ImIl L (s) over 
the invariant mass s of the final-state hadrons as follows: 



* - ia *f S t 1 - w) K 1 + 2 s?) imnT(s) + imnl(s) 

The appropriate combinations of two-point correlation functions are 



U J (s) = \V ud \' K dy (s)+U J , A (s) 



(2.5) 



(2.6) 



with \ Vij \ the corresponding Cabibbo-Kobayashi-Maskawa (CKM) quark mixing fac- 
tors. 

Experimentally, one can disentangle vector from axial-vector Dirac structures 
in the Cabibbo-allowed (ud) flavour decays. Since G-parity is not a good quan- 
tum number for modes with strange particles, this separation is problematic for the 
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Cabibbo-suppressed decays. It is then convenient to decompose the predictions for 
R T into the following three categories: 



R T = R t V + R t A + R t S . (2.7) 

Where R T y and R T) a correspond to the the two terms proportional to |V^| 2 in (2.6) 
and R T) s contains the remaining \V US \ 2 contributions. 

Exploiting the analytic properties of IT J (s), we can express (2.5) as a contour 
integral in the complex s-plane running counter-clockwise around the circle \s\ — M 2 : 

* - - L* t i 1 - w) 3 { 3 1 1 + flL+T(s) +iDL(s) ] ■ (2 - 8) 

We have used integration by parts to rewrite R T in terms of the logarithmic deriva- 
tives of the relevant correlators, 

D L+T {s) = -s A [n^( S )] , D\s) = ^ A [ s U L (s)} , (2.9) 

which satisfy homogeneous renormalization group (RG) equations. In this way, one 
eliminates unwanted (renormalization-scheme and scale dependent) subtraction con- 
stants, which do not contribute to any physical observable. 

For large enough —s, the contributions to D J (s) can be organized using the OPE 
in a series of local gauge-invariant scalar operators of increasing dimension D = 2n, 
times the appropriate inverse powers of — s. This expansion is expected to be well 
behaved along the complex contour \s\ — M 2 , except in the crossing point with the 
positive real axis [11]. As shown in eq. (2.8), the region near the physical cut is 
strongly suppressed by a zero of order three at s = M 2 . Therefore, the uncertainties 
associated with the use of the OPE near the time-like axis are very small. Inserting 
this series in (2.8) and evaluating the contour integral, one can rewrite R T as an 
expansion in inverse powers of M 2 [8], 



R T 



\Vud\ 2 + \V US \ 2 



s 



EW 



\l + 5' EW +5^+ J2 (cos 2 Be +sin 2 9 C 

[ D=2,4,- 



(2.10) 



where sin 2 9c = |14s| 2 /[|^4d| 2 + |Kts| 2 ] an d we have pulled out the electroweak cor- 
rections Sew = 1-0194 [12] and 5' EW ~ 0.0010 [13]. 

The dimension-zero contribution 5^ is the purely perturbative correction, ne- 
glecting quark masses, which, owing to chiral symmetry, is identical for the vector 
and axial-vector parts. The symbols 8\f^ = [5\fy + <^a]/2 stand for the average 
of the vector and axial-vector contributions from dimension D > 2 operators; they 
contain an implicit suppression factor 1/M®. 
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A general analysis of the relevant S^ V , A contributions was presented in ref. [8]. 
A more detailed study of the perturbative piece was later performed in ref. [14], 
where a resummation of higher-order corrections induced by running effects along 
the integration contour was achieved with RG techniques. More recently, the leading 
quark-mass corrections of dimension two have been investigated in ref. [6]; these 
contributions are the dominant SU(3) breaking effect, which generates the wanted 
sensitivity to the strange quark mass. 

In order to simplify the presentation, we will relegate a detailed compilation of 
the different contributions to R T to the Appendix. 



3. Moments of the Hadronic Invariant Mass Distribution 

The measurement of the invariant mass distribution of the final hadrons provides 
additional information on the QCD dynamics. The moments [15] 

can be calculated theoretically in the same way as R T = R®°. The corresponding 
contour integrals can be written as 

R k T l = -m<f — hFt l +T {x) D L+T {M 2 T x) + 4F?{x) D L {M 2 T x)\ , (3.2) 

where all kinematical factors have been absorbed into the kernels 

T kl +T {x) = 2 (1 - xf +k 

y y g U n (Q + k + n) + 2(3 + k + n)x 

h(l-n)\n\ { ' (3 + k + n)(4 + k + n) ' ^ 



rt\x) = 3(1- xf +k V - \, , } , ' . (3.4) 
LV ' V ; ^ (l-n)\n\ 3 + k + n K J 

Table 1 shows the explicit form of these kernels for the moments which we are going 
to analyze in the following sections. 

One can rewrite the moments R kl as an expansion in inverse powers of M 2 , 
analogously to (2.10). The corresponding contributions from dimension D operators 
will be denoted by 8 k f D \ 



4. SU(3) Breaking 

The separate measurement of the Cabibbo-allowed and Cabibbo-suppressed decay 
widths of the r [4] allows one to pin down the SU(3) breaking effect induced by the 
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(k,l) 






(0 0) 


(1 - x) 3 (1 + x) 


(1 - x) 3 


(1,0) 


£ (i _ x f (7 + 8x ) 


I (I--) 4 


(2,0) 


JL (i _ X )5 (4 + 53.) 


§(! -xf 


(1,1) 


±(l-a;) 4 (l + 2a;) 2 


!(l-x) 4 (l + 4a;) 


(1,2) 


2io (1 - a;) 4 (13 + 52a; + 130a; 2 + 120a; 3 ) 


^ (1 - a;) 4 (1 + 4x + 10x 2 ) 



Table 1: Explicit values of the relevant kinematical kernels. 



strange quark mass, through the differences 

6b E W 2^ [°ud - d us 



T>kl 

§R kl = K,V + A 



T)kl 
K r,S 



D>2 



(4.1) 



These observables vanish in the SU(3) limit, which helps to reduce many theoretical 
uncertainties. In particular they are free of possible (flavour-independent) instanton 
and/or renormalon contributions which could mimic dimension two corrections. 

4.1 Dimension— Two Corrections 

The dimension-two corrections to the r hadronic width are perturbative contri- 
butions proportional to m 2 . We studied in a previous paper [6] the uncertainties 
associated with these corrections, in the limit m u = m d = 0. The main conclusion of 
that work was the relatively large uncertainty in the prediction of these corrections 
arising from the very bad behaviour of the J = L component. We give here the 
general result, for arbitrary light quark masses. 

In terms of the running quark masses and the QCD coupling, the D = 2 contri- 
butions to the correlation functions take the form 



D=2 



An 2 .i 



rn. 



'i-es)+m*(-es)} Y.c L n +T (0a n (-es) 



n=0 



±m, t (-es) mji-es) £ et +T (0 a^s) 

71=1 

E/n +T (0«"(-^)l, (4.2) 
n=2 ) 



E^(-£ 2 *) 

L k 



Df jy/A {s) 



D=2 



8tt 2 M2 



mi {-es)Tm 3 {-es)] 2 £ d L n (t)a n (-es), (4.3) 



n=0 



where a = a s /7i and £ is an arbitrary scale factor of order unity. The coefficients 
c^ +T (£), e^ +T (£), fn +T {i)i an( l ^n(0 are constrained by the RG equations satisfied 
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by the corresponding D J (s) functions (in fact all of them follow the same RG equa- 
tions). Their scale dependence is given in Appendix B, together with the values of 
the presently known coefficients in the MS scheme. 

Inserting these expressions into the contour integral (3.2), the corresponding 
D = 2 contributions to Rf, 5®yJ A , are given by analogous expansions, with the 
running coupling a n (— £ 2 s) replaced by the functions: 



f>kl(n) ( v -1 / dx T ki . 



x) 



m(-£ 2 M 2 x) 
m(M 2 ) 



-£ 2 M 2 *) , (4.4) 



B 



kl (n) . 



dx 



2m J\x\=i x Ly 1 



m( 



-i 2 M 2 T x 



1 2 



a n {-eM 2 T x) . 



(4.5) 



m(M 2 ) 

These integrals only depend on = a s (£ 2 M 2 )/7r, log(£), and the expansion 
coefficients $ and 7j of the QCD beta and gamma functions. They were already 
studied in ref. [6] for the case (k,l) = (0,0). 

We only need the contribution to 5R^ 1 , which is given by 



SR. 



m 2 (M 2 ) , 2 x (2) 
M£>ew — r— 2 — \L - e d ) l\ kl [a T ) , 



D=2 M 2 

where a T = « s (M 2 )/tt, e d = m d /m s = 0.053 ± 0.002 [16] and 1 



(4.6) 



a m K) = j 3 £ £ L n +T (0Bl^(as) + £ d L n (0Bl tw (a,) 

4 I n=0 n=0 J 

= \{ZA L k + T (a T ) + A L kl (a T )} . (4.7) 

The longitudinal series A kl (a T ) is unfortunately quite problematic. The bad 
perturbative behaviour od Df jV/A {s) gets reinforced by running effects along 

J ' 1 D=2 

the integration contour, giving rise to a badly defined series. The convergence can 
be improved [6] by fully keeping the known four-loop information on the function 
integrals Bj l ^ n \a^), i.e. using in eqs. (4.4) and (4.5) the exact solution for m(— £ 2 s) 
and a(— £ 2 s) obtained from the RG equations. This "contour-improved" prescription 
[14] allows us to resum the most important higher-order corrections, but the resulting 
"improved" series is still rather badly behaved. For instance, 

A&)(0.1) = 1.5891 + 1.1733 + 1.1214 + 1.2489 + • • • (4.8) 

which has 0(a 2 ) and 0(a 3 ) contributions of the same size. On the contrary, the 
J = L + T series converges very well: 

A^ +T (0.1) = 0.7824 + 0.2239 + 0.0831 - 0.0000601 c% +T +■■■ (4.9) 



1 Notice that is slightly different from the analogous quantity A L+T defined in ref. [6] 

with d£ +T (0 = S£ +T (£) + fn +T (0- Thc SU(3) singlet component f£ +T (0 drops out in 5R k T l . 
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Fortunately, the longitudinal contribution to A kl (a T ) is parametrically sup- 
pressed by a factor 1/3. Thus, the combined final expansion looks still acceptable 
for the first few terms: 

Ag(O.l) = 0.9840 + 0.4613 + 0.3427+ (0.3122 - 0.000045 4 +T ) +••• (4.10) 

Nevertheless, after the third term the series appears to be dominated by the longi- 
tudinal contribution, and the bad perturbative behaviour becomes again manifest. 
Using cf +T ~ C2 +T (c2 +T /cf +T ) ss 323, the fourth term becomes 0.298; i.e. a 5% 
reduction only. We can then take the size of the 0(a 3 ) contribution to A kl as an 
educated estimate of the perturbative uncertainty. 

The final numerical values of the relevant perturbative expansions are shown in 
Table 2. We have used the value of the strong coupling constant determined by the 
total hadronic decay width [3]: 

a s (M?) = 0.35 + 0.02. (4.11) 

Two different errors are quoted in Table 2. The first one gives the estimated theo- 
retical uncertainties for the central value of a s (M^), while the second one shows the 
changes induced by the present uncertainty in the strong coupling. 

Since the longitudinal series (4.8) seems to reach an asymptotic behaviour at 
0(a 3 ), we have taken the following criteria in our numerical estimates. The central 
values of A$ (a T ) have been evaluated adding to the fully known 0(a 2 ) result one half 
of the 0(a 3 ) contribution. The 0(a 3 ) running effects in the L+T contribution have 
been also included; the remaining 0(a 3 ) contribution from the unknown constant 
C3 +T was estimated above to be less than 1% in Aqq . To estimate the associated 
theoretical uncertainties, we have taken one half of the size of the last known per- 
turbative contribution plus the variation induced by a change of the renormalization 
scale in the range £ G [0.75, 2] (added in quadrature). Finally the central values have 
been obtained by symmetrizing the error bars. 



(M) 




A£,(a T ) 




(0,0) 


0.97 + 0.10 + 0.07 


5.1 + 2.1 + 0.5 


2.0 + 0.5 + 0.1 


(1,0) 


1.37 + 0.12 + 0.05 


5.3 + 2.5 + 0.7 


2.4 + 0.7 + 0.1 


(2,0) 


1.70 + 0.30 + 0.09 


5.8 + 3.2 + 0.8 


2.7+1.0 + 0.2 


(1,1) 


-0.37 + 0.11 + 0.02 


-0.45 + 0.66 + 0.25 


-0.39 + 0.25 + 0.07 


(1,2) 


0.02 + 0.03 + 0.01 


0.23 + 0.17 + 0.08 


0.07 + 0.06 + 0.02 



Table 2: Numerical values of the relevant D = 2 perturbative expansions for a s (M 2 ) = 
0.35 ± 0.02. The first error shows the estimated theoretical uncertainties taking a s = 0.35; 
the second one shows the changes induced by the present uncertainty in the strong coupling. 
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Notice from Table 2 that the L + T series is more sensitive to the value of the 
moment parameter k than the L series. On the other side, the two last moments with 
/ 7^ give rise to perturbative expansions for A^(a T ) which are clearly unreliable; 
therefore, we will discard these two moments in our final m s analysis. 

4.2 Dimension-Four Corrections 

The SU(3)-breaking piece of the D = 4 contribution to the correlation functions is 
given by 



s 2 



D L u $ +a {s)-D l u $ +a {s)] d _ a = -A50 4 (-es) £ «)« n K S «) 

n=0 

+ (1-3) E{(l + 3) ^ +T (0-3^ +T (0}«"(-^),(4.12) 

71 71=0 1 ' ' 



sM 2 



K d y +A (s) - Dt y+ As)\ D=A = 2S0 4 (-es) 



-- 2 mt(-es) (1-4) E {(1 + 3) p£(0 +Jn(0] (4-13) 

^ n=0 k 

[2^(0- 3^(0 +m)]}a n (-es), 

where e u = m u /m s = 0.029 ± 0.003 [16], was defined before, and 

50 4 (/i 2 ) = (0|m s ss-m d J<i|0)(/i 2 ). (4.14) 

The quark condensates are defined in the MS-scheme, at the scale fj 2 = — £ 2 s. 
The perturbative expansion coefficients g£ +T (0, ^ +T (0> <?n +T (0; 
and J^(£) are given in Appendix C. 

Inserting these expressions into the contour formula (3.2), one gets the corre- 
sponding contributions to 8R}f. They can be written in the form 



5R k T l 



d _ 4 = 12S ew {3^r^(l -3) [{l + e 2 d )T kl {a T )-2elS kl {a T ) 

-^ T ^QuM). (4-15) 

The normalization of the perturbative expansions Tki(a T ), Sm{cl t ), and Qki{oj T ) has 
been chosen so that, for the lowest-order moments, these quantities are just one at 
leading order, i.e. T o(0) = £00 (0) = Qoo(0) = 1. Their explicit expressions are given 
in Appendix E, for the (k, I) values which are going to be relevant in our analysis. 
Table 3 here and Tables 9, and 10 in Appendix E show their corresponding numerical 
values. 

In principle, the SU(3)-breaking condensate SO^M 2 ) could be extracted from 
the r decay data, together with m S) through a combined fit of different SRf moments. 
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(k I) 




Si An, l 


OiAa ) 

V kl 


(0 0) 


1.5 ± 0.5 ± 0.1 


1.0 ± 0.2 ± 0.05 


1.08 ± 0.03 ± 0.01 


(1,0) 


3.2 ±0.7 ±0.2 


1.9 ±0.6 ±0.2 


1.52 ±0.03 ±0.01 


(2,0) 


5.2 ±0.9 ±0.3 


3.5 ± 1.4 ±0.2 


1.93 ±0.02 ±0.003 


(1,1) 


-2.0 ±0.3 ±0.2 


-1.5 ±0.3 ±0.1 


-0.41 ±0.02 ±0.005 


(1,2) 


0.45 ±0.20 ±0.20 


0.3 ±0.1 ±0.1 


-0.02 ±0.01 ±0.003 



Table 3: Numerical values of the relevant D = 4 perturbative expansions for a s (M 2 ) = 
0.35 ± 0.02. The first error shows the estimated theoretical uncertainties taking a s = 0.35; 
the second one shows the changes induced by the present uncertainty in the strong coupling. 



However, this is not possible with the actual experimental accuracy. In the future 
this could be the best determination of the SU(3)-breaking condensate <50 4 (M 2 ). 

We can estimate the value of <50 4 (M 2 ) using the constraints provided by the 
chiral symmetry of QCD. To lowest order in Chiral Perturbation Theory [17], £0 4 (/x 2 ) 
is scale independent and is fully predicted in terms of the pion decay constant and 
the pion and kaon masses: 



5o 4 (m t 2 ; 



= (m s - m d ) (0|gg|0) ~ -jl (m 2 K - m 2 ) ~ -1.9 x 10~ 3 GeV 4 

0( P 2) ' " V 1 



(4.16) 



Here, (0|gg|0) denotes the quark condensate in the chiral limit, which we take to 
be approximately given by 2m(0|gg|0) ~ m (0\dd + uu\0) ~ —j 2 m\ [17, 18], where 
rh = (m u ± m d )/2. 

We can improve this estimate, taking into account the leading 0(p 4 ) corrections 
through the ratio of quark vacuum condensates 2 , 

• <4 - i7) 

This ratio has been phenomenologically estimated to be around 0.6 ~ 0.9 for scales 
between 1 and 2 GeV where the scale dependence is very soft [19, 20, 21]. To 
be conservative, we have enlarged slightly its allowed range to include the SU(3) 
symmetric value v s — 1. Taking this correction into account, we get our final estimate 



50 4 (M T 2 



= [v s m s -m d ) (0\dd\0) ~ -— i (v s - e d ) / 2 m 2 
o( p 4) 2m 

~ -(1.5 ±0.4) x 10~ 3 GeV 4 , (4.18) 



2 Strictly speaking this ratio is UV scale dependent in QCD. As shown in Appendix D, this 
dependence is canceled by m 4 terms and is then of order p 8 in the momentum expansion. For typical 
hadronic scales, these 0(m 4 ) corrections are very small and will be neglected in the following. 
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where we have used the known values of quark mass ratios [16] m s /rh = 24.4 ± 1.5 
and td given before. At this order, SO^fi 2 ) is still scale independent in QCD. 

For typical hadronic scales the quark condensate gives a sizeable contribution 
to SRf, proportional to -4tt 2 50 4 (M 2 )/M 4 m (5.9 ± 1.6) x 1(T 3 , which is much 
larger than the remaining 0(m 4 ) corrections. Those m 4 contributions are of the 
same order than the scale dependence of <50 4 (/i 2 ), which is smaller than the accuracy 
of our estimate (4.18). To be consistent, we will therefore neglect all m 4 corrections 
together with the scale dependence of SO^fj, 2 ). 



4.3 Higher Dimension Corrections 

The leading order coefficients of dimension six and eight corrections for two point 
functions have been studied in the MS scheme [22, 23, 24, 25, 26, 27, 28, 29]. However, 
in view of the theoretical uncertainties in the dimension two and four corrections and 
the unknown values of the dimension six and eight condensates, we shall not include 
the D > 6 contributions and regard them as an additional theoretical uncertainty. 

To get an order of magnitude estimate of the possible size of those effects, let us 
neglect their logarithmic dependences and parameterize the leading D = 6 contribu- 
tions to the two-point correlators as 



s 

s 2 M 2 



D=6 
D=6 



-3<50 6 L+T , 



2 50%. 



Dtd,v+A( s ) ~ D^ s V+A (s) 
The corresponding contribution to the different 5R* 1 moments is then: 

12tt 2 



(4.19) 



5R° T ° 



5R1 



5R 20 



5R? 
5R 12 



D=6 



D=6 



D=6 



D=6 



D=6 



M r 6 
12tt 2 
M T 6 
12tt 2 
M T 6 

127T 2 

— 50 L+T 
MS 6 • 



3 50% +T -4 50% 



3 50% +T 



[2 50% +1 



6 50% 



8 50% 



(4.20) 



[50% +T + 2 50, 



One could expect that the leading D = 6 contributions come from four-quark 
operators, because they are not suppressed by light quark masses (the G 3 operator 
is flavour symmetric and therefore does not contribute to 5Rf). Obviously, only the 
L + T piece gets such a contribution, which can be obtained from ref. [8] in the 
vacuum saturation approximation, 



5R 



oo 



-S 



D=6 



EW 



36tt 2 



50^ T » S 



EW 



256vr 4 (0|ss|0) 2 (M 2 ) - (0|Jd|0) 2 (M 2 ) 



M r 6 



10 



o c X 00 ( D = 6 ) ~~ V s 

6 b EW d ud — - — 



(0.6 ±2.3) x 10" 



In that case, we also have 



5 Ft 



.00 



SR 



10 



D=6 



D=6 



35R 



ii 



D=6 



-36 R. 



12 



D=6 



D=6 



(4.21) 



(4.22) 



To get the final number, we have used the measured value of the Cabibbo-allowed 
correction [30, 31] S, 



00 (D=6) 
ud 



0.001 ±0.004. 



The size of dimension six corrections proportional to four-quark operators is 
smaller than the uncertainty from the dimension four corrections for the three mo- 
ments that we are going to use [(k,l) = (0,0), (1,0), (2,0)]. We thus conclude that 
dimension six and higher corrections are negligible given the actual experimental 
accuracy (see Table 4) and do not add any further uncertainty to the m s determina- 
tion at present if one uses the three moments above. They will become eventually 
important for the determination of the SU(3)-breaking condensate <50 4 (M^). 

Notice that higher-order corrections could be important and even dominant in 
some cases. For instance, in the moment (k,l) = (1,2) there is a strong suppression 
of the contributions with dimensions two and four (see Tables 2 and 3), which makes 
necessary to consider the D = 6 terms. 



5. Numerical Analysis 



Discarding 0(m 4 ) corrections, which are much smaller than the present experimental 
uncertainties, and up to dimension six corrections, which are around eight times 
smaller than the uncertainty in the dimension four contribution for the three moments 
[(k, I) = (0, 0), (1, 0), (2, 0)] considered here, 



m 2 s (M 2 T ) 



Ml 



2(l- e i)AgV) [12S EW 



5R: 



kl 



+ 4tt 2 \\ a T ' Q k i{a T ) 



M r 4 



(5.1) 



The ALEPH collaboration has measured [4] the weighted differences 8Rf for five 
different values of (k,l). The experimental results are shown in Table 4, together 
with the corresponding m s (M 2 ) values. Since the QCD counterparts to the moments 
(k,l) = (1,1) and (1,2) have theoretical uncertainties larger than 100 %, we only use 
the moments (k, I) = (0,0), (1,0), and (2,0). 

The experimental errors quoted in Table 4 do not include the present uncertainty 
V us \. To estimate the corresponding error in m s , we take the following numbers 



m 



published by ALEPH [4, 30]: R% +A = 3.486 ±0.015, R° T ° S = 0.1610 ±0.0066, \V X 



0.9751 ±0.0004 and IK 



0.2218±0.0016. This gives 3 SR™ 



ud 



0.394±0.135±0.047, 



3 Using the PDG98 [1] values (with the constraint \V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1) \V ud \ = 0.97525± 
0.00046 and \V US \ = 0.2211 ± 0.0020, we get SR™ = 0.372 ± 0.136 ± 0.060, which lowers 5 MeV the 



central value of m s {M 2 ) from 5R°°. 
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(k I) 


5R kl 


m (M 2 ) (MeV) 


(o 0) 


0.394 ± 0.137 


143 ± 31 ± 18 


(1,0) 


0.383 ±0.078 


121 ± 17 ± 18 


(2,0) 


0.373 ± 0.054 


106 ± 12 ±21 


(1,1) 


0.010 ±0.029 




(1,2) 


0.006 ±0.015 





Table 4: Measured [4] values of 5R!f and the strange quark mass at M T . The first error 
is experimental and the second is from the uncertainty in the QCD counterpart. 

where the second error comes from the uncertainty in \V US \ an d translates into an 
additional uncertainty of 10 MeV in the strange quark mass. Since the ALEPH 
collaboration does not quote the values of R* l v+A and R* l s for the other moments, 
we will put the same \V US \ uncertainty to the other two moments. 

Taking the information from the three moments into account, we get our final 
result for m s (M 2 ): 

m s (M 2 ) = (119 ± 12 ± 18 ± 10) MeV = (119 ± 24) MeV (5.2) 

The first error is experimental, the second reflects the QCD uncertainty and the 
third one is from the present uncertainty in \V US \. Since the three moments are 
highly correlated, we have taken the smaller individual errors as errors of the final 
average. 

At the usual scales /i = 1 GeV and = 2 GeV (used as reference values by 
the QCD Sum Rules and lattice communities, respectively), our determination (5.2) 
corresponds to 

m s (l GeV 2 ) = (164 ± 17 ± 25 ± 14) MeV = (164 ± 33) MeV (5.3) 

and 

m s (4 GeV 2 ) = (114 ± 12 ± 17 ± 10) MeV = (114 ± 23) MeV . (5.4) 

6. Phenomenological Subtraction of the J = L Piece 

In order to avoid the large theoretical uncertainties associated with the bad pertur- 
bative behaviour of A^(a T ), it would be nice to have "experimental values" for the 
J = L + T contributions to 5R}f . 

Using the positivity of the spectral functions, the known pion and kaon poles 
provide the lower bounds 

bn 14(a) > 27r/ 2 5( S -m 2 ), lmIl L us {s) > 2vr f K 5(s - m\) , (6.1) 
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(k I) 


5R k K , rr. 


mJM 2 ) (MeV) 




< 0.411 


< 287 


(1,0) 


< 0.351 


< 246 


(2,0) 


< 0.326 


< 202 


(1,1) 


< 0.030 




(1,2) 


< 0.020 





Table 5: Bounds on 5R kl L+T and the strange quark mass at M T . 



which translate into upper limits on the corresponding J = L contributions to R^. , 

2+k / \ 1+2 



R k T ] L = -247T 



M r dS 

M 2 
dx 



Ml 



Ml 



Imn L (s) 



= _ 4vn I ™ jf ( x )D L (M 2 x). 

J\x\=l X ' 



(6.2) 



After subtracting the Goldstone boson pole, Im Hfj(s) is proportional to light quark 
masses squared. Since m s ^> m u ^, we can then safely conclude 



8Rf L > 48tt 2 



Pk 



X 



m 



2 \ 
K \ 



m 



i+i 



K 



Ml 



Ml 



Ml 



J TT 

Ml 



o \ 2+k / o \ 

Ml ) \Ml 



. (6.3) 



Subtracting this contribution from the measured values of 5R k , one gets upper 
bounds on 5R kl L+T . Therefore, using the relation 



IM 2 
ml(M 2 )- ~ 



1 



3(1 - e 2 ) A L k + T (a T ) [ 12S EW 



$ R t\l+T 



+ 47T : 



,S0 4 (M 2 ) 



M* 



Q 



L+T, 
kl \ 



(6.4) 



we can get non-trivial upper bounds on m s (M 2 ). The resulting numerical values are 
given in Table 5. 

To improve on these bounds it would be necessary to have a better understanding 
of the J = L spectral functions. 



7. Comparison with the ALEPH Analysis 

The ALEPH collaboration has also performed a phenomenological analysis of the 
weighted differences SRf in Table 4. However, the resulting values of the strange 
quark mass quoted by ALEPH [4] are larger: 

m (M' 2 ) — I (M9 ± 6 fit) MeV (inclusive) , 

' ;i76^ 8 7 :-l^±8 fit ±ll J=0 )MeV (L + Tonly). 
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To derive these numbers, ALEPH has used our published results in references [6], 
[8], and [14]. Since we have analysed here the same ALEPH data with improved 
theoretical input, it is worthwhile to understand qualitatively the main origin of the 
numerical difference. 

ALEPH makes a global fit to the five measured moments. Moreover, they 
also fit two additional parameters Sf^ and trying to extract higher-order non- 
perturbative corrections from the data. As we have pointed out before, the last two 
moments have theoretical errors (in the leading D = 2 contribution) larger than 
100% and therefore are unreliable. Unfortunately, the sensitivity to the small D > 6 
corrections comes precisely from these two moments (see our comments at the end 
of Section 4.3), which makes the fitted Sf' 8 ^ values very questionable. 

Owing to the asymptotic behaviour of A^(a T ), ALEPH truncates the contribu- 
tion of this perturbative series at 0(a) (in the "inclusive" analysis), neglecting the 
known (and positive) 0(a 2 ) and 0(a 3 ) longitudinal contributions. But the relevant 
quantity for this analysis is A^(a r ) which does not reach its miminum term before 
0(a 3 ) [6] [see Eq. (4.10)], so it is inconsistent to neglect the large and positive longi- 
tudinal contributions to A^f (a r ) as ALEPH did. Thus, they use a smaller value of 
A^i(a T ) and, therefore, get a larger result for the strange quark mass, because the 
sensitivity to this parameter is through the product m 2 s (M^) A^(a T ). Since they put 
rather conservative errors, their result [the value quoted as "inclusive" in eq. (7.1)] 
is consistent with ours. Nevertheless, it is a clear overestimate of m s because 4 they 
underestimate A kl (a T ). 

In order to avoid the large perturbative corrections in the longitudinal piece, the 
ALEPH collaboration has made a second analysis, subtracting the J = L contri- 
bution in a way completely analogous to the one presented in Section 6. However, 
besides subtracting the pion and kaon poles, ALEPH makes a tiny ad-hoc correction 
to account for the remaining unknown J = L contribution, and quotes the resulting 
number as a m s (M 2 ) determination [the value quoted as "L + T only" in eq. (7.1)]. 
Since they add a generous uncertainty, their number does not disagree with ours. 
It is clear, however, from our discussion in Section 6, that this is actually an upper 
bound on m s (M 2 ) and not a determination of this parameter. 

8. Summary 

We have analysed the SU(3) breaking effects in the semi-inclusive r hadronic width 
in complete generality. This has been used to obtain the strange quark mass from 
the recent ALEPH measurement of the inclusive Cabibbo-suppressed decay width 

4 In fact, larger values of m s were also obtained in the first determinations of this parameter from 
QCD sum rules, performed at 0(a s ). Once the large higher-order perturbative corrections to the 
corresponding correlators were known, the resulting m s values shifted down by a sizeable amount. 
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and several moments of its invariant mass distribution. We get 

m s (M 2 ) = (119 ± 12 ± 18 ± 10) MeV = (119 ± 24) MeV (8.1) 

to 0(a 3 s ) within the MS scheme 5 . The first error comes from the experimental un- 
certainty, the second one from the uncertainty in the QCD counterpart and the third 
from the uncertainty in \V US \. 

At the customary scales where quark masses are quoted, this result translates 

into 

m s {\ GeV 2 ) = (164 ± 17 ± 25 ± 14) MeV = (164 ± 33) MeV (8.2) 

and 

m s (4 GeV 2 ) = (114 ± 12 ± 17 ± 10) MeV = (114 ± 23) MeV . (8.3) 

This agrees within errors with the findings in ref. [6], where only 5R®° was used. 

Subtracting the known kaon and pion poles, we have also obtained an upper 
bound on the strange quark mass, 

m s (M 2 ) < 202 MeV, (8.4) 

which corresponds to m s (lGeV 2 ) < 277 MeV and m s (4GeV 2 ) < 194 MeV. This 
bound is completely free from the problems associated with the bad perturbative 
behaviour of the J = L contribution. Our result is compatible with the lower bounds 
presented in refs. [32, 33, 34]. 

There is a great deal of activity calculating the strange quark mass by the lattice 
community [35]. The results are still very confusing and the spread in values obtained 
using different approximations to QCD is quite large. For a critical view of the 
situation see [36]. 

The latest QCD Sum Rules determinations [37, 38, 39, 40, 41, 42] have been 
obtaining results which are very compatible with our number. The systematic error 
in those determinations is however still unclear. 

The sum of the up and down quark masses has been determined with Finite 
Energy Sum Rules in ref. [18] with the result (m u + m d )(4GeV 2 ) = (9.8 ± 1.9) 
MeV. Using the ratio of light quark masses 2m s /{m u + m d ) = 24.4 ± 1.5, obtained 
within 0(p 4 ) Chiral Perturbation Theory and large N c [16], this result is also in nice 
agreement with our determination 

We have not made any attempt to reduce the theoretical error, which we defer to 
a future publication. As stated in ref. [6] , once the invariant mass distribution of the 
final r-decay hadrons is known, it should be possible to find weighted distributions 

5 There is a remaining 0(a 3 ) contribution from the unknown constant Cg +T , which was estimated 
in Section (4.1) to modify m s by less than 0.5 %. 
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with smaller theoretical errors for the dimension-two QCD counterpart. If the SU(3)- 
breaking data from tau decays is improved at future facilities, it could be the source 
of precise determinations of both the strange quark mass and the SU(3)-breaking 
condensate 80 4 (M^). 
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A. Dimension— Zero Corrections 

Though they have been extensively studied in refs. [8] and [14], for the sake of 
completeness, we give here also the dimension-zero corrections to R* 1 . They are 
flavour independent and are identical for the vector and axial-vector correlators. 
Moreover, only the transverse piece gets D = contributions: 



m s Dtf /M d __ o Av 



l -Y,IUQa"(-es). (A.l) 



The coefficients K n (g) are constrained by the homogeneous RG equations satis- 
fied by the Adler function D(s): 



for n > 1, and 



Thus, 



£ A K n (0 = J2(k-n)[3 k K n ^ k (0 , (A.2) 
a ^ k=l 



±K (O=0. (A.3) 



K (0 = K , 
K 1 (£) = K 1 , 
K 2 (S) = K 2 -/3 1 K 1 hiS, 

K 3 (0 = K 3 - [f3 2 K x + 2^K 2 ] In f + 0{K X In 2 i , 

K 4 (0 = K 4 ---- (A.4) 
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The factors (3k are the expansion coefficients of the QCD beta function, defined 

as 

(3(a) = ^ = £/^ (A.5) 

which are known to four loops [43, 44]. For three flavours and in the MS scheme, 

9 3863 n 140599 445 

0i = — , 02 = -8 , 03 = , At = Cs ■ A -6 

2 ^ ' ^ 192 ' ' 2304 16 v ; 

The perturbative expansion of the Adler function is fully known up to order a z s . At 
£ = 1, its coefficients have the values [45, 46]: 

299 _ T ^ 58057 779 75 , 
= * = A- 2 = — -DCs, ^, = — -—01-^. (A.7) 

The perturbative component of R T is given by 

^ = EW4 n) W> (A.8) 

n=l 

where the functions [14] 

4?K) = -^rW«1-»), (A.9) 

are contour integrals in the complex plane which only depend on ag = a s (£ 2 M^)/7r, 
ln(£), and the expansion coefficients 0i of the QCD beta function. A detailed analysis 
of this contribution to R T and its associated uncertainty can be found in ref. [14] . 

B. Dimension— Two Corrections 

The dimension-two corrections to the correlators Df- v , A {s) are given in eqs. (4.2) 
and (4.3). The corresponding expansion coefficients obey the RG equation 

Z^CniO = £ [2 7fe — (n — k) k ] (7„_ fe (0 , (B.l) 
a ^ k=l 

for n > 1 and 

^G)(0=0, (B.2) 

where the generic notation (7(£) stands for c^ +T (£), e^ +T (£), fn +T (C)i an d 
Therefore, 

Cb(0 =Co, 



17 



as 



C 1 (0 = C 1 + 2 7l C lne, 

C 2 (0 = C 2 + [2 72 C + (2 7l - ln£ + 7i (27i - Pi)C In 2 £ , 

^3(0 = C 3 + [2 73 Co + (2 72 - / 3 2 )C 1 + 2( 7l - Pi)C 2 ] hi£ 

+ [(- 7l /3 2 + 2 72 (2 7l - /3 1 ))C + (71 - A)(27i - ln 2 £ 
+ ^ 7 i( 7 i- / 9i)(2 7l -A)C ln 3 e, 
(74(0 = C 4 + -- - • (B.3) 
The factors 7fc are the expansion coefficients of the QCD gamma function, defined 

which are known to four loops [44, 47, 48]. For three flavours and in the MS scheme, 

91 8885 r> 

7i = 2 , 72 = — , 73 = -^gg- - 5Cs , 

2977517 9295 135 125 
74 = ^0736" " W Cs + ~ C4 " "6" C5 - (R5) 

The J = L + T coefficients are known to 0(a 2 s ) [49, 50, 51, 52, 53, 54, 55] while 
the J = L coefficients are known to 0(a^) [49, 50, 51, 55]. Their values at £ = 1 are 

23077 179 520 

432 + ~54 ^ 3 ~ ~W^ 5 ' 
769 55 5 

~54 ~~ 27 ^ 3 ~ 2T 5 ' 
-| + 2C„ (B.6) 

and 

9631 35 

■ "144" ~ y ^ 3 ' 

In the limit m u = = 0, taken in ref. [6], the expansion of the J = L + T 
correlator is governed by the combinations 4 +T (£) = c^ +T (£) + fn +T {C)- 

C. Dimension— Four Corrections 

The dimension-four corrections to the correlators D^ v ^ A (s) can be written in the 
form: 

D ^v/a(s)\ d _ a = y 2 E ^ +T (^ s ) ° n (-^ s ) > (C- 1 ) 



c 


1, 


c i 


13 

~T ' 


c 2 _ 


e 


0, 


e i 


2 

~ 3' 


e 2 _ 


fL+T _ 

Jo 


0, 


fi +T 


= 0, 


rL+T _ 
J2 




4 


= 1, 


= 


3 ' 




4 


474£ 
51 


>953 
34 


91519 
216 ^ 3 
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D 




+ 7 ^(m tT m J ) 2 ]T n L n (Z,s)a n (-es)), (C.2) 



n=0 



where 





(m! + mj) h L n (0 ± fntrnj ~k L n {t) + Y, m l • (C4) 



In these expressions, the running masses and vacuum condensates are defined in the 
MS scheme at the scale ji 2 = — £ 2 s, i.e. 



m t = rmi-fs) , (G 2 ) ( = (0\G 2 \0) (-£ 2 s) , (mOM)* = (0\m t q jqj \0) (-£ 2 s) . 



In ref. [8] the D = 4 contributions were given in terms of scale-invariant con- 
densates. This simplifies the R T contour integration, but introduces inverse powers 
of a s in some m 4 terms, generating larger quark-mass corrections which cancel nu- 
merically with the condensate contributions [56]. With the minimally subtracted 
operators (C.5), used here, one gets slightly more stable numerical results for the 
dimension-four mass corrections. 

The quark condensate contribution to the longitudinal correlator (C.2) is fixed 
to all orders in perturbation theory by a Ward identity. 



The perturbative expansion coefficients in eqs. (C.l) and (C.2) are known to 
0(a 2 ) [27, 29, 57, 58, 59, 60, 61, 62, 63, 64] for the condensate contributions, 



(C.5) 




r L+T (0 = 0, 





1, P L 2 +T (0 = l 




r T (0 = o, 






(C.6) 



while the m 4 terms have been only computed to 0(a) [24, 62, 63, 65, 66]: 
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h^ +T (0 = 1 + In $ , ^ +T (0 = y - 2 C 3 + f In $ + 4 In 2 £ , 
S^ r (0 = l, ^ +T (0 = y-^3 + 8 1ne, 

j L+T (0 = o, jT +t (0 = o, 

^(0=0, 5^(0 = 0, (C.7) 

-41 28 
Jtf (0 = 1 + In £ , h{ (0 = - - 2 C 3 + y In t + 4 In 2 £ , 

# (0 = 1 + I In £ , fcf (0 = § - \ C 3 + y In £ + \ In 2 £ ■ 

io L (0 = o, jT(0 = o. (C.8) 

The scale dependence of all these coefficients is fixed by the homogeneous RG 
equations satisfied by the corresponding D J (s) functions: 

a ^ k=l a, 5 

a ^ fe=l J fc=l 

^^ T (0 = |7ipf +T (0, 

e^*n> 2 (0 = E(*-») = 0, (c.9) 



H _ n _ 1 n 

£y ^(0 = E [47, - (n - *) A] " ^5*) ^(0 , 

U, S fc=l Z fc=0 

4^ +T (0 = -^7So)r T (0, 

j n— 1 o n— 1 

£y fe T (0 = E [47, — (n — k) (3 k \ k^) - ± ^ {k) , 

a ^ k=l k=0 

^*f +T (o = ~75o ) ^ +T (o, 

^ 5& T (0 = E [47, - (n - *) A] ^J(0 , ^ ^ L+T (0 = 
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J n— 2 i n— 1 

fei£ftO = E[47,-(n-A;) / 3 fc ]j^(0 + -E^owfe T (0 

^ — 1 /c — 1 

~ 7 E 7o(fc) r n-fc (0 > 77 ^n<l (0 = 0) 

^ fc=0 a, 5 

4 # +T (0 = ^o(d ?r T (0 - ^§ 0) rt +T (o , 

J n— 2 i rt— 1 

a, 5 fc=l z ^fc=2 

^^(0 = 0, (cio) 

H ~ 1 n 

£77 ^>i(0 = -77§n) + E - (n- k) (3 k \ h L n _ k {0 , 

a, 5 z fc=l 

£ 7^ (0 = — 2 7o(o) > 

d ~ 1 n 

^ ^>i(0 = -3 7o ( n) + E [47, — (n — k) (3 k \ k L n _ k {t) , 

^ ' k=l 

£ 7^ *o(0 = _ g 7o(o) > 

d ~ n ~ 2 
£77 #>s(0 = E [47, - (n - k) ft] -jtAt) , 
^ fe=i 

^5»< 2 (0 = o. (en) 

D. Scale Evolution of the D = 4 Operators 

The factors 7^ appearing in the RG equations (C.10) and (C.ll) are the anomalous 
dimensions of the QCD vacuum energy, 

4 E = -&<;+ E rn k q k q k , (D.l) 

k=u,d,s 

where 9^ is the trace of the QCD energy-momentum tensor, in the three light flavour 
effective theory, and 

dE 3 

^"rf = 7^ E m i m ) 7o J ( a ) ! 7o J (a) = 7o(o) + a 7o(i) + ' • ' ( D - 2 ) 

" i,j=u,d,s 

The anomalous dimension matrix 70(a) is symmetric, i.e. Jo (a) = 7o*( a )- More- 
over, since QCD is flavour blind, 7" (a) = 7^ (a) and all non-diagonal elements Y(f j (a) 
are equal. To two loops, Jo (a) is proportional to the identity with [45] 

7§o) = -2, 7S(i) = -| » ( D - 3 ) 
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for i = u,d, s. The first non-diagonal terms appear at three-loops. 

The anomalous dimension matrix 70(a) governs the scale evolution of the D = 4 
operators. After using the QCD equations of motion, there are three types of gauge- 
invariant operators of dimension four, namely, G 2 = Gf^Gfy, rriqq, and m 4 . In 
minimal subtraction-like schemes, these three operators mix under renormalization 
[45]: 



% «?> = -. 



d(3(a) 



(G 2 } + — — a £ in-, "I 

i,j=u,d,s 



da v ' 4tt 2 . 1 3 da 



2 2 dYo{a) 97(a) „ _ 



da 



=u,d,s 



(D.4) 



(D.5) 



u— m 4 = — 4 7(a) m 4 . 



(D.6) 



One can introduce [45] the following scale-invariant condensates, which are com- 
binations of the previous minimally subtracted operators: 



^(aG^j = (3(a) (G^G$) + 47(a) £ (m« 

i=u,d,s 



— ]T m 2 m 2 7o », 

i,j=u,d,s 



(D.7) 



{miqjqj)i = (mqfij) 



3 mi rrij 



To(o) 



+ 



7o(i) 7o(o) #2 + 472 



47i 4 7 i /3i + 471 



47r 2 a + 471 

3m;mf f 53 1 _ , 

/,From these two invariants, one inmediately gets the scale evolution of the cor- 
responding MS operators. In particular, the SU(3)-breaking difference of quark 
condensates (4.14) satisfies: 



*0 4 ( M 2 ) = ,50 4 (M T 2 ) + -5j (l - £ J) mJ(M r 2 ) f ' 



a T + 



(D.9) 



m 4 (/i 2 ^ 



1 



1 — — a(a 2 ) 
24 VP y 



m 4 (M 2 ) a(/i 2 ) 
E. D = 4 Expansion Coefficients for 

The functions Q kl {a) = Q^(a) + Q^ +T (a), T kl {a) = T&(a) + T^ +T {a) and S kl (a) = 
Ski( a ) + Ski +T ( a ) i n ec f (4-15) are given by the following contour integrals 



Qki 



6m J\x\=i x 



rkl 50 A {-eM 2 x) 



dx 



o0 4 (M 2 ) 
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To be consistent with our estimate of 50 A(fJ. 2 ), which is scale independent, we will 
also neglect the <50 4 (/z 2 ) scale dependence in the functions Q J kl (a). In the following 
tables we summarize the known values of the different D = 4 perturbative expansions 
for the (k, I) values used in the paper. We also give their numerical values obtained 
form the exact integration along the closed contour, using the corresponding RG 
equations. 



(k,l) 


T kl (a) 


S k i(a) 


Qki(a) 


(0,0) 
(1,0) 
(2,0) 

(1,1) 
(1,2) 


l + (f-4C3-H« 

§[l + (f -fCa-Ha] 
f [l+df-ICa-l^a] 

I [l + m 


1 + fa 

2[l + (f -§Cs)a] 

-1 [1 + (¥ " Ks) «] 
±a 


2 1 — 3 a + 72 a 
2 1 — 2° 96~ a2 

-^a 2 

160 



Table 6: Known terms of the relevant D = 4 perturbative expansions. 
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(1,0) 
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/ 496 
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24 



15" 

- ^7T 2 
17" 
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_1 

4 
J_ 
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1 + fa 
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1 + 



1 
1 



fc/V 



37 
3 

41 
3 



1 + fa 



Q 



1 
1 
1 





Table 7: Known terms of the relevant D = 4 perturbative expansions for J = L. 
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(M) 



(0,0) 
(1,0) 
(2,0) 

(1,1) 
(1,2) 



1 + ffl 



l + (^f + 6C3 + 7T 2 ) a 



13 
12 

_IZ 
12 



-i , / 266 48 a 

+ (# 

40 



13 ^ 13 

24 a 4 
S3 — 



7r 2 ] a 



17 



17 



7T 2 I a 



1 + l^a 
i -t- 1Q5 a 
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Table 8: Known terms of the relevant D = 4 perturbative expansions for J = L + T. 



(M) 


TfczK) 


Skl( a r) 




(0,0) 


3.10 ±0.40 ±0.04 


2.1 ±0.7 ±0.2 


1.0 ±0.0 ±0.0 


(1,0) 


3.80 ±0.40 ±0.08 


2.7 ± 1.0 ±0.2 


1.0 ±0.0 ±0.0 


(2,0) 


4.4 ±0.3 ±0.2 


3.2 ± 1.3 ±0.2 


1.0 ±0.0 ±0.0 


(1,1) 


-0.66 ±0.07 ±0.06 


-0.6 ±0.2 ±0.1 


0.0 ±0.0 ±0.0 


(1,2) 


-0.03 ±0.10 ±0.07 


-0.01 ±0.02 ±0.01 


0.0 ±0.0 ±0.0 



Table 9: Numerical values of the relevant D = 4 perturbative expansions for a s (M 2 ) = 
0.35 ± 0.02 for J = L. The first error shows the estimated theoretical uncertainties taking 
a s = 0.35; the second one shows the changes induced by the present uncertainty in the 
strong coupling. 



(M) 




s L k ria T ) 




(0,0) 


-1.60 ±0.06 ±0.1 


-1.1 ±0.5 ±0.1 


0.08 ±0.03 ±0.01 


(1,0) 


-0.6 ±0.4 ±0.2 


-0.6 ±0.4 ±0.1 


0.52 ±0.03 ±0.01 


(2,0) 


0.80 ±0.60 ±0.05 


0.25 ±0.05 ±0.05 


0.93 ±0.02 ±0.003 


(1,1) 


-1.30 ±0.20 ±0.15 


-0.86 ±0.18 ±0.02 


-0.41 ±0.02 ±0.005 


(1,2) 


0.48 ±0.20 ±0.20 


0.3 ±0.1 ±0.1 


-0.02 ±0.01 ±0.003 



Table 10: Numerical values of the relevant D = 4 perturbative expansions for a s (M 2 ) = 
0.35 ± 0.02 for J = L + T . The first error shows the estimated theoretical uncertainties 
taking a s = 0.35; the second one shows the changes induced by the present uncertainty in 
the strong coupling. 



References 

[1] Review of Particle Physics, C. Caso et al., Eur. Phys. J. C 3 (1998) 1 

[2] Proceedings of the Fifth Workshop on Tau Lepton Physics -TAU 98- (Santander, 
14-17 September 1998), A. Pich and A. Ruiz (eds), Nucl. Phys. B (Proc. Suppl.) 76 
(1999). 



24 



[3] A. Pich, Tau Physics, in "Heavy Flavours II", A.J. Buras and M. Lindner (eds), 
Advanced Series on Directions in High Energy Physics - Vol. 15 (World Scientific, 
Singapore, 1998), p. 453 [hep-ph/9704453] 



[5] 

[e; 

[7] 
[8] 
[9 

[io; 
[ii 

[12 

[13 
[14 
[15 
[16 



[17 

[is; 

[19 

[20; 

[21 



R. Barate et al. (ALEPH), hep-ex/9903015; hep-ex/9903014 

A. Hocker, Nucl. Phys. B (Proc. Suppl.) 55C (1997) 379; M. Davier, Nucl. Phys. B 
(Proc. Suppl.) 55C (1997) 395; S. Chen, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 256; 
S. Chen, M. Davier, and A. Hocker, Nucl. Phys. B (Proc. Suppl.) 76 (1999) 369 

A. Pich and J. Prades, J. High Energy Phys. 06 (1998) 013; J. Prades and A. Pich, 
Nucl. Phys. B (Proc. Suppl.) 74 (1999) 309; J. Prades, Nucl. Phys. B (Proc. Suppl.) 
76 (1999) 341 

K.G. Chetyrkin, J.H. Kiihn, and A. A. Pivovarov, Nucl. Phys. B 533 (1998) 473 

E. Braaten, S. Narison, and A. Pich, Nucl. Phys. B 373 (1992) 581 

S. Narison and A. Pich, Phys. Lett. B 211 (1988) 183 

E. Braaten, Phys. Rev. Lett. 60 (1988) 1606; Phys. Rev. D 39 (1989) 1458 

E.C. Poggio, H.R. Quinn, and S. Weinberg, Phys. Rev. D 13 (1976) 1958 

W. Marciano and A. Sirlin, Phys. Rev. Lett. 61 (1988) 1815; Phys. Rev. Lett. 56 
(1986) 22 

E. Braaten and C.S. Li, Phys. Rev. D 42 (1990) 3888 

F. Le Diberder and A. Pich, Phys. Lett. B 286 (1992) 147 
F. Le Diberder and A. Pich, Phys. Lett. B 289 (1992) 165 

H. Leutwlyer, in Proc. of the NATO Advanced Study Institute on Masses of Fun- 
damental Particles (Cargese, 1996), M. Levy, J. Illiopoulos, R. Gastmans, and J.-M. 
Gerard (eds), NATO ASI Series B:Physics, Vol. 363 (Plenum Press, New York, 1997), 
p. 149 [hep-ph/9609467]. 

J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465; Ann. Phys. (NY) 158 
(1984) 142 

J. Bijnens, J. Prades, and E. de Rafael, Phys. Lett. B 348 (1995) 226, J. Prades, 
Nucl. Phys. B (Proc. Suppl.) 64 (1998) 253 

S. Narison, QCD Spectral Sum Rules, World Scientific Lecture Notes in Physics , Vol. 
26 (1989). 

H.G. Dosch, M. Jamin, and S. Narison, Phys. Lett. B 220 (1989) 175 
S. Narison, Phys. Lett. B 358 (1995) 113 



25 



[22] L.-E. Adam and K.G. Chetyrkin,P/i?/s. Lett. B 329 (1994) 129 

[23] L.V. Lanin, V.P. Spiridonov, and K.G. Chetyrkin, Sov. J. Nucl. Phys. 44 (1986) 892 

[24] M. Jamin and M. Miinz, Z. Physik C 60 (1993) 569 

[25] M. Jamin and M. Miinz, Z. Physik C 66 (1995) 633 

[26] D.J. Broadhurst and S.C. Generalis, Phys. Lett. B 165 (1985) 175 

[27] S.C. Generalis and D.J. Broadhurst, Phys. Lett. B 139 (1984) 85 

[28] W. Hubschmid and S. Mallik, Nucl. Phys. B 207 (1982) 29 

[29] S.C. Generalis, J. Phys. G 16 (1990) 367 

[30] ALEPH Coll., R. Barate et al., Eur. Phys. J. C 4 (1998) 409 

[31] OPAL Coll., K. Ackerstaff et al., Eur. Phys. J. C 7 (1999) 571 

[32] L. Lellouch, E. de Rafael, and J. Taron, Phys. Lett. B 414 (1997) 195 
L. Lellouch, Nucl. Phys. B (Proc. Suppl.) 63 (1998) 272 
E. de Rafael, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 258 

[33] F.J. Yndurain, Nucl. Phys. B 517 (1998) 324 

[34] H.G. Dosch and S. Narison, Phys. Lett. B 417 (1998) 173 

[35] D. Becirevic, Report on light quark masses. Plenary talk at Lattice' 99 (July 1999, 
Pisa, Italy) 

R.D. Kenway, Nucl. Phys. B (Proc. Suppl.) 73 (1999) 16 
[36] C.R. Allton, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 242 

[37] P. Colangelo, F. De Fazio, G. Nardulli, and N. Paver, Phys. Lett. B 408 (1997) 340 
[38] M. Jamin, Nucl. Phys. B (Proc. Suppl.) 64 (1998) 250 

[39] C.A. Dommguez, L. Pirovano, and K. Schilcher, Phys. Lett. B 425 (1998) 193 

[40] K. Maltman, Phys. Lett. B 428 (1998) 179; Phys. Rev. D 58 (1998) 093015 

[41] K. Maltman, The Strange Quark Mass from Finite Energy Sum Rules, 
hep-ph/9904370 

[42] S. Narison On the Strange Quark Mass from e + e~ Data and Test of the SU(2) Isospin 
Symmetry, hep-ph/9905264 

[43] T. van Ritbergen, J.A.M. Vermaseren, and S.A. Larin, Phys. Lett. B 400 (1997) 379 

[44] K.G. Chetyrkin, B.A. Kniehl, and M. Steinhauser, Phys. Rev. Lett. 79 (1997) 2184 

[45] V.P. Spiridonov and K.G. Chetyrkin, Sov. J. Nucl. Phys. 47 (1988) 522 



26 



[46] S.G. Gorishny, A.L. Kataev, and S.A. Larin, Phys. Lett. B 259 (1991) 144 

[47] J.A.M. Vermaseren, S.A. Larin, and T. van Ritbergen, Phys. Lett. B 405 (1997) 327 

[48] K.G. Chetyrkin, Phys. Lett. B 404 (1997) 161 

[49] L.R. Surguladze, Sov. J. Nucl. Phys. 50 (1989) 372 

[50] S.G. Gosrishny, A.L. Kataev, S.A. Larin, and L.R. Surguladze, Mod. Phys. Lett. A 5 

(1990) 2703 

[51] S.G. Gosrishny, A.L. Kataev, S.A. Larin, and L.R. Surguladze, Phys. Rev. D 43 

(1991) 1633 

[52] S.G. Gorishny, A.L. Kataev, and S.L. Larin, Nuovo Cim. 92 (1986) 119 

[53] W. Bernreuther and W. Wetzel, Z. Physik C 11 (1981) 113 

[54] K.G. Chetyrkin and J.H. Kiihn, Phys. Lett. B 406 (1997) 102 

[55] K.G. Chetyrkin, Phys. Lett. B 390 (1997) 309 

[56] K.G. Chetyrkin and A. Kwiatkowski, Z. Physik C 59 (1993) 525 

[57] M.A. Shifman, A.L Vainshtein, and V.I. Zakharov, Nucl. Phys. B 147 (1979) 385, 
448 

[58] G.T. Loladze, L.R. Surguladze, and F.V. Tkachov, Phys. Lett. B 162 (1985) 263 

[59] C. Becchi, S. Narison, E. de Rafael, and F.J. Yndurain, Z. Physik C 8 (1981) 335 

[60] L.R. Surguladze and F.V. Tkachov, Nucl. Phys. B 331 (1990) 35 

[61] E. Bagan, J.I. Latorre, and P. Pascual, Z. Physik C 32 (1986) 43 

[62] K.G. Chetyrkin, S.G. Gosrishny, and V.P. Spiridonov, Phys. Lett. B 160 (1985) 149 

[63] S.C. Generalis, J. Phys. G 15 (1989) L225. 

[64] P. Pascual and E. de Rafael, Z. Physik C 12 (1982) 127 

[65] S.C. Generalis, J. Phys. G 16 (1990) 785. 

[66] D.J. Broadhurst and S.C. Generalis, Improved Two-Loop Mass Corrections: The Role 
of the Quark Condensate, Open University preprint OUT-4102-22 (1987) (unpub- 
lished) . 



27 



